Electromagnetic response of high-T c superconductors — 
the slave-boson and doped-carrier theories 
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We evaluate the doping dependence of the quasiparticle current and low temperature superfluid 
density in two slave-particle theories of the tt't" J model - the slave-boson theory and doped-carrier 
theory. In the slave-boson theory, the nodal quasiparticle current renormalization factor a vanishes 
proportionally to the zero temperature superfluid density ps(0); however, we find that away from 
the ps(0) — > limit a displays a much weaker doping dependence than ps(0). A similar conclu- 
sion applies to the doped-carrier theory, which differentiates the nodal and antinodal regions of 
momentum space. Due to its momentum space anisotropy, the doped-carrier theory enhances the 
value of a in the hole doped regime, bringing it to quantitative agreement with experiments, and 
reproduces the asymmetry between hole and electron doped cuprate superconductors. Finally, we 
use the doped-carrier theory to predict a specific experimental signature of local staggered spin cor- 
relations in doped Mott insulator superconductors which, we propose, should be observed in STM 
measurements of underdoped high-T c compounds. This experimental signature distinguishes the 
doped-carrier theory from other candidate mean-field theories of high-T c superconductors, like the 
slave-boson theory and the conventional BCS theory. 



I. INTRODUCTION 



Motivation 



The phenomenology of high-temperature supercon- 
ducting (SC) cuprates is most striking when electron 
occupancy per unit cell is close to unity 1 . In this 
underdoped regime, the SC critical temperature (T c ) 
vanishes as the interaction-driven Mott insulator is ap- 
proached despite the strong binding of electrons into 
Cooper pairs 2,3 . This deviation from BCS theory is re- 
flected in the anomalous metallic pseudogap state, whose 
Fermi surface appears to be partially gapped 4,5 . The 
physics that determines the value of T c in underdoped 
cuprates is a highly debated 6-18 fundamental question 
which lacks a fully satisfactory answer. 

In this context, it is important to sort the specific 
roles played by phase fluctuations of the order param- 
eter and by fermionic quasiparticle excitations in de- 
stroying superconductivity at finite temperatures. Since 
the zero temperature superfluid phase stiffness ps(0) de- 
creases together with T c 19 ~ 24 , phase fluctuations are cer- 
tainly detrimental to long-range phase coherence 6 , as ev- 
idenced by the magnetic vortices observed in the normal 
state 25 . Yet, to reconcile T c with the measured bare 
phase stiffness 26 requires that thermally excited nodal 
quasiparticles reduce ps, thus promoting vortex prolifer- 
ation at a temperature lower than expected from phase- 
only arguments 7,10,12 . The important role played by the 
d x 2_ y 2-weLve quasiparticles receives further experimental 
support from penetration depth measurements in under- 
doped YBCO films which display characteristic quasi- 
particle behavior, namely the linear in T suppression of 
ps, up to T > T c /2 (sec, for instance, Fig. 1 of Ref. 
20). Remarkably, this T-linear regime extends all the 



way to T c in severely underdoped samples which, in ad- 
dition, violate the T c oc ps(0) relation applicable in pure 
phase-fluctuation models 21-24 . The above strongly sup- 
ports that the T c scale in underdoped cuprates is set by 
the effective parameters ps(ty and dps/dT. Simultane- 
ously describing the above two parameters in consistency 
with experiments is the main problem we address in this 
paper. 

In the remaining of this introductory section, we il- 
lustrate that present theoretical understanding cannot 
easily reconcile the experimentally observed behavior of 
p s {0) with that of dp s /dT (Sec. IB). The latter pa- 
rameter reflects the interaction induced renormalization 
of the quasiparticle current and, as such, in Sec. I C 
we discuss how interactions are expected to renormalize 
the quasiparticle current throughout the entire Brillouin 
zone. We also note that these renormalization effects 
can be probed by scanning tunneling microscopy (STM) 
experiments, which can provide valuable information to 
distinguish between different quasiparticle descriptions of 
underdoped curpates. Finally, in Sec. ID we present the 
layout of the full paper. 



B. The effective parameters ps(0) and dps/dT 

The effective parameter ps(0) quantifies the coupling 
between an applied gauge field and the superconducting 
condensate, as evidenced in the Meissner effect. The de- 
pletion of ps (0) as the Mott insulator is approached 19-24 
is in stark contrast with the prediction from BCS the- 
ory, namely ps(0) oc 1 — x, where x is the density of 
carriers doped away from half-filling. This sharp de- 
viation follows from the large Coulomb repulsion that 
suppresses (enhances) charge (phase) fluctuations and, 
indeed, the observed ps(0) ~ x behavior is captured 
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by the slave-boson (SB) theory of the tJ model 12 ' 27,28 , 
which is a microscopic approach that explicitly imple- 
ments the suppression of charge fluctuations. The same 
behavior is encountered in theories of phase fluctuating 
d-wave superconductors 11,14 ' 15,29 ' 30 , whose effective the- 
ory is similar to that of the SB approach 11 ' 14 . However, 
the above success in reproducing ps{0) ~ x is not ac- 
companied by a similar fate when it comes to addressing 
dps/dT's experimental results. In fact, the microscopic 
SB approach predicts too small values of dps / dT in the 
limit x — > 0. 

The parameter dps/dT is often characterized in 
terms of the nodal current renormalization factor a = 
[(2-kva/vf In 2) dps/dT] 1 / 2 , where vf and va arc the 
Fermi and gap velocities respectively. In the BCS theory, 
a = 1. In the SB theory, however, the effect of Coulomb 
repulsion leads to a ~ Ps(0) in the limit x — > o 7 ' 9 ' 11,12 . 
This is commonly regarded as a major setback since ex- 
periments show that a vanishes sublincarly in ps(0) 22 ~ 24 . 
In addition, it casts doubt on the applicability of the SB 
formalism to simultaneously describe how SC quasipar- 
ticles and the SC condensate couple to an applied elec- 
tromagnetic field. 

We remark that the above mismatch between the SB 
theory and experiments occurs in the limit x — > 0, in 
which case the SB theory ignores the emergence of the 
antiferromagnetic (AF) phase. Therefore, in this paper 
we extend previous work in the literature that uses the 
slave-particle framework to address the supcrfluid density 
ps in the limit x — > 7 ~ 9 , and calculate the low energy 
and long wavelength electromagnetic response function 
of a doped Mott insulator superconductor away from the 
above limit. We specifically consider two slave- particle 
theories of the tt't" J model, namely the SB and the 
doped-carrier (DC) theories 27,28 ' 31,32 , for which we calcu- 
late the nodal current renormalization factor a as a func- 
tion of x. We argue that, in this respect, slave-particle 
theories may compare to experiments better than often 
thought. In fact, we find that both the SB and the DC 
slave-particle approaches predict that, for x > 0.05 and 
in the considered parameter range 2 J < t < 5 J, the dop- 
ing dependence of a is much weaker than that of ps(0), 
in agreement with underdoped cuprates' data 20 . 



C. Renormalized quasiparticle current distribution 

As we state above, superfluid density measurements 
probe the quasiparticle current renormalization at the 
nodal points of d x 2_ y 2-w&ve superconductors. Interac- 
tions, however, also renormalize the current of quasipar- 
ticles away from the nodes, an effect which should mani- 
fest itself in experiments, as we overview in what follows. 

We know that a finite supercurrent J s shifts the su- 
perconducting quasiparticle dispersion E k and, to linear 
order in J s , we have 



where we introduce the vector potential A to represent 
the supercurrent J s — psA. [Note that in Eq. (1) we set 
the speed of light to c = 1; in what follows, we also take 
the electric charge to be e = 1, as well as h = I.] The 
(hole) quasiparticle current j k characterizes how excited 
quasiparticles affect the superfluid density ps and, in the 
BCS theory, it is given by the expression 



Jk = — 



de 



N 



dk 



(2) 



E k {A) = E k (0) -j k ■ A 



(1) 



where is the normal state (hole) energy dispersion 
which, in the present discussion, wc approximate by = 
2t[cos(k x ) + cos(ky)]. [Up to a constant scale factor, Fig. 
1(a) resembles this quasiparticle current distribution.] 

We note that, according to the BCS theory, the quasi- 
particle current j k is completely determined by the nor- 
mal state dispersion (this applies all the way deep into 
the superconducting state). Since, at many levels, the 
phenomenology of overdoped cuprates is compatible with 
the BCS theory, we still expect Eq. (2) to hold in these 
samples, with given by the appropriate free elec- 
tron dispersion. However, the free electron dispersion 
should not be used to calculate the quasiparticle cur- 
rent of superconducting underdoped cuprates. This then 
raises the question of what normal state dispersion we 
should use to calculate the quasiparticle current distri- 
bution j k in underdoped cuprates. The answer to this 
question may come from angle-resolved photoemission 
spectroscopy (ARPES) experiments in half-filled cuprate 
compounds showing that the single hole dispersion is 
roughly given by e^ F oc [cos^fe^) + cos(2k y )] 33 , in which 

case the quasiparticle current j k = — [see Fig. 1(f)] 
considerably differs from the BCS-like result in Fig. I (a). 
This suggests that, perhaps, in the cuprates' underdoped 
regime we should use a quasiparticle current derived from 
a dispersion that interpolates between ej^ and e k F . 

From the above discussion, we see that the quasipar- 
ticle current distribution j k is an important quantity 
that can reveal new characteristics of underdoped high- 
T c superconductors which lie beyond the BCS paradigm. 
Hence, it is relevant to calculate and predict j k using 
different approaches to the high-T c problem. It is also 
significant to experimentally measure j k , as such a mea- 
surement could prove to be instrumental in either ruling 
out or validating candidate theories to the high-T c prob- 
lem. 

As a step in this direction, below we calculate the 
quasiparticle current distribution using two different ap- 
proaches, namely the aforementioned SB and DC ap- 
proaches. We find that these yield quite distinct distri- 
butions of the quasiparticle current throughout the Bril- 
louin zone [see Figs. 1(a) and 1(c)]. The mean-field SB 
approach gives rise to a quasiparticle current distribution 
which, in the absence of intra-sublattice hopping pro- 
cesses, is essentially the BCS result multiplied by the cur- 
rent renormalization factor a. The DC approach results 
in a quasiparticle current distribution which, instead, in- 
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tcrpolates between the BCS and AF current distributions 

de N de AF 

3k = ~~dt and 3k = respectively. 

In this paper, we propose that measuring the tunnel- 
ing differential conductance from a metal tip into the 
superconducting C u O plane in the presence of a super- 
current provides a way to distinguish the above quasi- 
particle current distributions. Specifically, we calculate 
how the tunneling differential conductance is affected by 
an applied supercurrent in the C u O plane, and find that 
the different quasiparticle current distributions lead to 
different supercurrent dependences of the tunneling dif- 
ferential conductance (see Fig. 4). We further argue that 
this effect may be probed by STM experiments, which 
thus could distinguish the DC theory description of high- 
T c superconducors from alternative theories (particularly 
those that ignore the momentum space differentiation of 
the nodal and antinodal regions, such as the SB theory 
and the conventional BCS theory). 

D. Paper layout 

The paper is organized as follows. In Sec. II we in- 
troduce the formalism used to calculate the low energy 
and long wavelength electromagnetic response function 
within the SB and DC frameworks. In order to ob- 
tain the non-zero temperature electromagnetic response 
function in the static and uniform limit we follow Ioffe 
and Larkin and resort to the random-phase approxima- 
tion (RPA), which accounts for the effect of fluctuations 
around the mean-field saddle point up to Gaussian level. 
This calculation yields the doping dependent quasiparti- 
cle current and ps{T) low temperature behavior in the 
SB and DC theories, whose results we discuss and com- 
pare in Sec. III. In this section, we also study the para- 
metric dependence of ps(T) on t/J values. In the DC 
framework we use in this paper, the momentum space 
anisotropy follows from the role of high-energy AF corre- 
lations between local moments in superconducting doped 
Mott insulators. Hence, the aforementioned comparison 
between the SB and the DC theory results identifies the 
effect of short-range AF correlations in the electromag- 
netic response of a <i-wave superconductor close to a Mott 
insulator transition. We conclude that AF correlations 
enhance (suppress) the nodal quasiparticle current in the 
hole (electron) doped regime of cuprate superconductors. 
In the hole doped case, and for t/J = 3, that enhance- 
ment leads to 0.5 < a < 0.6 when 0.1 < x < 0.2, in quan- 
titative agreement with experiments 35-37 . In this dop- 
ing range, the DC theory quasiparticle-driven T c scale 
T^ p = p s {0)/(dps/dT) ~ J/10 is an order of magnitude 
lower than in the SB theory and, in addition, it agrees 
with the cuprates' T c scale. Our results also show that 
the intriguing weak temperature dependence of ps{T) in 
electron doped compounds 38,39 is consistent with the ob- 
served momentum space anisotropy 40 , which we consider 
to follow from the strong local AF correlations. Since the 
formalism introduced in Sec. II also allows one to study 



the coupling between the SC quasiparticles and an ap- 
plied supercurrent, in Sec. IV we use this fact to predict 
a specific experimental signature of a traversing super- 
current in the tunneling differential conductance that, 
we propose, should be detected in STM measurements. 



II. FORMALISM 

We set to calculate the low energy and long wavelength 
electromagnetic response function of a doped Mott in- 
sulator superconductor. In particular, we consider two 
different families of slave-particle wave functions, one de- 
scribed by the SB d-wave SC ansatz 27,28 and the other by 
the DC (i-wave SC ansatz 31 ' 32 . In addition, we choose the 
energetics to be given by the tt't" J model Hamiltonian 

H t j = E J v S *- S 3 - E *y( 5 I,<A<7 + h.c.) (3) 

(ij {ij),° 

where c\ a = c\ a (l — c\ _ er Cj ) _ (T ) are the Gutzwiller pro- 
jected electron operators, Si — c\crci are the electron 
spin operators and <x are the Pauli matrices. Also, 
Jij = J for NN sites, and tij — t,t',t" for first, sec- 
ond and third NN sites, respectively. Since the SC nodal 
quasiparticles are the sole gapless excitations, below we 
resort to the SB (Sec. II A) and DC (Sec. II B) mean-field 
theories. In both cases, we include the effect of gapped 
collective modes at the RPA level. 



A. Slave-boson framework 

We first determine how the well studied 7 ' 8 ' 12,27 ' 28 ' 41 SB 
d-wave superconductor couples to an applied electromag- 
netic field. In the SU(2) SB framework the projected 
electron operators are decoupled as — -j=ipjhi and 

cjj_ = -^i>T {—io~2)hi, where fa are the chargeless and 
spin-I/2 spinon fermionic operators in the Nambu rep- 
resentation, and hi are the spinless and charge-e holon 
bosonic operators 27 . If one rewrites Eq. (3) in terms of tpi 
and hi, and further applies the Hartree-Fock-Bogoliubov 
decoupling scheme, one obtains the quadratic SB mean- 
field Hamiltonian Hff F = H^ B + H^ B , where 42 : 

H " B = E ^f Tr + «o- (E ^) " 

^ 1 (4) 

(ij) 

Hf* = E y PnVji] +«o.(E h\*hi) 

{l3) t 1 (5) 

-W.E h l hi - E f (X^A + h.c) 

i (ij) 
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The holon chemical potential fih controls the density 
{h\hi) = x and the Lagrange multiplier a = a (T 3 im- 
plements the SU{2) projection constraint at mean-field 
level (ipjcripi + h\<jhi) = 0. In the d-wave SC ansatz 
h\ = [h 0], where h is the holon condensate magnitude, 



xcr 3 , U ii+& = x&3 + Acti, and U, 



ii+y = X&3 - A(Ti. 

Since only holons carry electric charge, minimally cou- 
pling an electromagnetic gauge field to Hf B F amounts 
to replacing h\Uijhj in Eq. (5) by h\Uij exp(iAij)hj. 
In order to account for the effect of fluctuations around 
mean-field we also must replace Uij and Vij in Eqs. 
(4) and (5) by [\Uij exp(ia,ij .cr) + i exp(iay .<r)Uij] and 



\\Vij exp(ia,ij.(r) + | exp(iaij.<r)Vij], respectively, where 
ttij is the SU(2) gauge field that describes collective 
modes in the SB framework 12,28,41 . In what follows, 
we consider the resulting minimally coupled Hamiltonian 
HfP F (Aij,a,ij) in the static and uniform limit, and thus 
recast A i+i a = A^u.x + A y u.y and a l+ m = a$u.x + 
a y u.y. 

If we ignore the contribution from collective modes, 
the electromagnetic current and response function 

dF SB (A,a) 



are 



jSB 
J A,u 



OA,, 



and 



ll AA,uv 



N dA a 9A S 



(A,a) 



A,a=0 

, where F SB (A, a) is the free-energy 

A,a=0 

obtained from Hf B F (Aij,a,ij). However, as shown by 
Ioffe and Larkin 34 , the above modes are important to cor- 
rectly determine how strongly correlated superconduc- 
tors couple to the electromagnetic field. In the SC state 
these modes are gapped and we only keep free-energy 
terms up to quadratic order in a fi . Integrating out a&, 
we obtain the electromagnetic current and response func- 
tion within RPA, namely, 



J 

n 



SB 



SB 



SB 



n 



SB 



(nff) 
nff) 



SB 
AA 



rSB 



yrSB 



(6) 
(7) 



where J^ B 



u sb 

Aa.uv 



dF SB (A,a) 



d 2 F 



da-c 
SB (A,a) 



A,a=0 



(u SB V 

\ 11 aA,vu J 



dA a dan 



A,a=0 



and II 



SB 



N 



d z F bB (A,a) 



<9a fl dat. 



A,a=0 



We note that and jf f cor- 

respond to the sums over momentum space Jf 2 ? = 

bE^wtgi and J ^ = wE b , k n b , k j^,l where 
rib^k is the occupation number of the mean-field sin- 
gle particle state in band b and with momentum fc, 
whose energy dispersion ef B {A, a) implies the quasi- 

particle current components j^'^? = 



dAu 

Following the above formula, 



and 



■ a,SB = deZ%(A,a) 
Jb.k.it 0a fl 

introduce the SB theory quasiparticlc current (within 
RPA) 

■SB _ 9Jj B _ .A,SB _ -f-rSB (ttSB}- 1 -a, SB 
Jb.k.u — _ Jb,k,u L1 Aa • \ ll aa ) ■ Jbte.u 



dnb,k 



(8) 



where we use J^ B = 0, which applies in thermal equi- 
librium and in the absence of external fields. Eq. (8) 



is equivalent to = -^jiil, where E^(A) 



e b,k (A a )la=-(n SB ) _ i n SB A 1S tne sm gl e particle disper- 
sion obtained in the presence of an applied gauge field 
Au after integrating out a&. Wc additionally obtain the 
SB theory supcrfiuid density (within RPA) p| B from 



n 



SB 



FS °u 



43 



B. Doped-carrier framework 

We now determine how an applied electromagnetic 
field couples to the DC c?-wave SC state in the static 
and uniform limit. The only difference when compared 
to the SB approach sketched in Sec. II A has to do 
with the specific decoupling of the Gutzwiller projected 
electron operators which, in the DC framework, read 32 

Here, s„ = 



(i + Sa sf) 4 

(+!),(-!) for a =T,|. di, a 



l 

V5 



r(l-d\ 



T di- a ) is the 

charge-e and spin-1/2 projected doped carrier operator 
(di, which has the same quantum numbers as the holes 
doped into the Mott insulator, has been called the do- 
pon operator in Rcf. 31). Further writing the above 
spin operators in terms of chargeless spin-1/2 spinons as 
Si = hfjcrfi leads to the DC mean-field Hamiltonian 



ttDC _ ttDC 
n MF — n i> 



H° C + H B g, where 31,32 : 



H * "16 



(ij)eNN i 



(ij)eNN 



(9) 



i v=2,3 u£u NN 



i i 

_ ^E E t» (4+itBioA + h.c) 



(10) 



(11) 



i,v u£uNN 



Above, ipi and r\i are the spinon and dopon operators 
in the Nambu representation, u = ±x, ±y, u = ±x ± y 
and u = ±2x,±2y for v = 1,2,3 respectively, /i^ is 
the chemical potential that sets the doped carrier density 
x. J = (1 - x) 2 J, h = t, t 2 = J + (1 - x/0.3)t' and 
<3 = J/2+(l— x/0.3)t", where J, t, t' and t" parameterize 
the tt't" J model Hamiltonian. In the <i-wave SC ansatz, 
Ua+x = X&3 + Acti, Ua+y = xcr 3 - A(Ti, B l0 = -b <7 3 
and Bn = —610-3. Furthermore, the Lagrange multiplier 
do = ao<J3 sets (ipjcripi) = 0. As shown in Rcf. 32, the 
SB and the DC formulations are related to each other 
since the local singlet state of a dopon di and a spinon 
fi in the DC approach corresponds to the holon in the 
SB approach. Consequently, the DC theory mean-field 
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^0 = (fidi) is equivalent to the SB holon condensate 
magnitude h. 

In this paper, we extend previous work on the DC 
framework and introduce the electromagnetic gauge field 
Aij , as well as the fluctuations around the Uij mean-field 
described by the SU(2) gauge modes (analogous to 
those in the SB formulation). Before constructing the 
minimally coupled Hamiltonian H® [F (A^ , ) , we sum- 
marize how the fields in Eqs. (9), (10) and (11) trans- 
form under U(l) electromagnetic gauge transformations 
and under SU (2) gauge transformations associated with 
the aij modes: (i) tpi is an on-site field which carries 
no electric charge and which is in the SU(2) fundamen- 
tal representation; (ii) rji is an on-site field which carries 
electric charge and is invariant to St/ (2) gauge transfor- 
mations; (iii) is an on-site field which carries electric 
charge and which is in the SU(2) fundamental represen- 
tation; (iv) Uij is defined on the lattice bonds, carries no 
electric charge and is in the SU(2) adjoint representation. 
As a result, H^ F (Aij, a^) is obtained from H^p by re- 



placing: (i) (^Ui- 
^<T 3 )exp(ia i: j.cr) - 



) in Eq. (9) by WU K 



±exp (i ai j.cr)(^-Ui. 



1 1 
2 



^3); (h) 



Vi+S^Vi in E q- (10) by r]l +s cxp (iA i+Si (T 3 )a 3 i]i; (iii) 

ri\ +i t B i0^i in E q- ( n ) b y 4+u ex V (iA l+Sl c7 3 )B M ip i . Sim- 
ilarly to Sec. II A, we further recast Ai + m = A±u.x + 
Ay u.y and Oj+sj = u.x + a y u.y. 

Given H^ F (Aij, ay), we can determine the DC mean- 
field free-energy F DC (A, a), from which the DC electro- 
magnetic current and response function follow. Using the 
aforementioned Gaussian approximation, we have 



jDC 
R DC 



7 DC 

ttDC 
VV AA 



nSf. (nSP)" 1 . j™ 



jjDC 



l aA 



(12) 
(13) 



where 

OF 



jDC _ _ j_ 
J Am ~ N 



dF"°(A,a) 



J, 



J_ 

N 



(n DC V = 

l *-*-aA,vu I 



A.a=0 



ttDC 
11 AA,uv 



A. o=0 
j_ d 2 F DC 
N dA a dA t 



DC 



(A,a) 



jtDC 

*- 3 -Aa,iiv 



J_ 

N 



a z F" L (A,q) 



dA a da t 



A,a=0 



A,a=0 

and 



n 



DC _ 



_1_ d z F L 
N 



A,a=0 



As in the SB ap- 



jDC _ 
J A/a — 



l_ ^ -A, DC 

2-~ib,k n b,kJb,k,ii, 



N 



and 



proach, we recast 
jDC - J_ v nuu i a ' DC 

U a,u — N l^b,k llb < k Jb.k,u ' 

gous to that in Sec. II A. The DC theory quasiparticle 
current (within RPA) then is 



where the notation is analo- 



aDC = 
Jb,k,u — 



dA, a 



(14) 



e?£ (A,a)\ a= _( Tl nc)-i n nc A is the 



where E^(A) - , bk 
quasiparticle energy renormalized by a^'s Gaussian fluc- 
tuations, and e^(A, a) is H^ F (Aij, ajj)'s eigenenergy 
for the quasiparticle state in band b and with momen- 
tum k (once again, we assume thermal equilibrium and 
the absence of external fields, in which case J®? = 0). 
Finally, the DC theory superfluid density (within RPA) 

Ps 



Pg C follows from 
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III. RESULTS 

In this section, we discuss our results for the electro- 
magnetic quasiparticle current and response function of 
superconductors described by the SB and DC formalisms. 
The main difference between these two frameworks is 
that the above mean-field DC approach captures the ef- 
fect of high energy and short-range staggered local mo- 
ment correlations in the low energy SC properties 31 ' 44 . 
Hence, below, we compare the SB and DC results to learn 
how local AF correlations renormalize the electromag- 
netic response of doped Mott insulator superconductors 
in the static and uniform limit. Naturally, the SC electro- 
magnetic response is determined by the underlying mean- 
field order parameters that define the SC phase, namely, 
/i, A ^ in the SB approach, and bo, A ^ in the DC 
framework. Therefore, in order to establish a meaning- 
ful comparison between the doping dependence of our 
SB and DC results, below we take the SB mean-fields h 
and A to be equal to the self-consistent DC mean-field 
parameters &o an d A, respectively. All other mean- field 
parameters are self-consistcntly determined within each 
approach. 



A. Quasiparticle current 

In Figs. 1(a) and 1(b), we show the momentum de- 
pendence of - in the x = 0.1 hole doped SC regime 
(henceforth defined by the model parameters tuD — 3J, 
t' HD = —2t" HD = — J) and in the x = 0.1 electron doped 
SC regime (henceforth defined by the model parameters 
tED = 3J, t' ED = —2t' FD = +J), respectively. Aside 
from a small magnitude difference, which is attributed 
to the larger holon condensation in the hole doped side, 
these current patterns are very similar to each other 
and are reminiscent of a quasiparticle dispersion approx- 
imately proportional to (cosfe^ +cosk y ). This result fol- 
lows from the spin-liquid correlations-driven renormal- 
ization of the intra-sublattice hopping parameters t' and 
i" 42 > 45 ; which only enter Hff F as part of the products xt' 
and xt" . 

The above SB results considerably differ, both in di- 
rection and magnitude, from those obtained in the DC 
approach [see Figs. 1(c) and 1(d) for the momentum 
dependence of j^u m tn c aforementioned parameter 
regimes]. For instance, in the hole doped regime, a vor- 
tex configuration appears in the jbk u vector map near 
(±7r/2, ±7r/2) [Fig. 1(c)] which is absent in the corre- 
sponding j^ku Pl°t [Fig- l( a )]- The DC quasiparticle 
current magnitude in this region of momentum space is 
also visibly larger. A similar vortex configuration (ab- 
sent in the SB theory) and enhanced quasiparticle cur- 
rent magnitude occur in the electron doped regime fi 
plot near (±tt, 0) and (0, ±tt) [Fig. 1(d)]. 

The clear anisotropy between the hole and electron 
doped regimes in the DC approach reflects the role of AF 
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FIG. 1: Quasiparticle current for: (a) x = 0.1 hole doped 
SC regime in SB approach; (b) x = 0.1 electron doped SC 
regime in SB approach; (c) x = 0.1 hole doped SC regime 
in DC approach; (d) x = 0.1 electron doped SC regime in 
DC approach; (e) x = 0.1, t = 3J and £' = t" = in DC 
approach; (f ) x — AF state. Current vectors in (a), (b), (c), 
(d) and (e) [but not (f)] are plotted with same vector scale. 



correlations, whose staggered pattern leaves the intra- 
sublattice hopping parameters t' and t" largely unrenor- 
malized. The effect of t' and t" in the presence of 
AF correlations has been extensively addressed in the 
literature 45-51 , which indicates that t' ~ — 2t" < lower 
the energy of AF correlations close to (±7r/2, ±7r/2), 
while t' ps — 2t" > have a similar effect close to (±7r, 0) 
and (0, ±tt). This trend supports that the differences 
between the above jff fi and j^u pl° ts are a manifes- 
tation of the underlying AF correlations captured within 
the DC approach. Indeed, the aforementioned vortices 
in Figs. 1(c) and 1(d) are clearly reminiscent of the 
x = AF state quasiparticle current plot in Fig. 1(f) 
[this current pattern is given by j£ F = —dk 1l E£ F , where 
we take E£ F oc [cos(2fc x ) + cos(2/c y )] 33 ]. Also, the large 
quasiparticle current close to (±7r/2, ±7r/2) in the hole 
doped regime, and close to (±7r,0), (0, ±7r) in the elec- 
tron doped regime, is also consistent with the well known 
enhancement of quasiparticle features in these regions of 
momentum space due to AF correlations 31,45,46,48-52 . Fi- 
nally, we note that, when t' = t" = 0, AF correlations do 
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FIG. 2: Doping dependence of (a) po, (b) &hd, and (c) T® p , 
for t/J = 2, 3, 4, 5 (see legend) and t , t" — t' HD ,t'l ID . Results 
are shown for both SB and DC theories, (d) Doping depen- 
dence of (cxed / olhd) 2 ■ Full [dashed] line in (a), (b), (c), and 
(d) plots t/J = 3 results in DC [SB] theory. 



not differentiate between the above two regions of mo- 
mentum space 45,53,54 , and the DC quasiparticle current 
pattern along the (0, 7r) — (tt, 0) resembles the SB results 
[Fig. 1(e)]. 



B. Super-fluid stiffness 

We now address the doping dependent behavior of 
ps(T) at low temperature, in both the SB and DC the- 
ories. In Figs. 2(a)-2(c) we take the hole doped regime 
parameters tf,t" = t' HD ,t'^ Dl and analyze the depen- 
dence on values of t/J that pertain to the physically rel- 
evant regime 2 < t/J < 5. Specifically, Fig. 2(a) depicts 
the scaled superfluid density po = g p (t/ J)ps(0), where 
the scaling function g p (t/J) = iJ/t in the SB theory and 
g p (t/J) = 1 in the DC theory. These results show that, in 
the considered parameter range, the SB approach yields 
ps(0) oc xt while in the DC approach ps(0) ~ xJ is al- 
most independent of t. This shows that AF correlations 
can renormalize the SC condensate kinetic energy scale 
from xt to x J. 

Fig. 2(b) depicts the scaled nodal quasiparticle cur- 
rent renormalization factor olhd = 9a(t/ J)olhd, where 
g a (t/J) = (SJ/t) 1 / 2 in the SB approach and g a (t/J) — 
Xt/ZJ) 1 / 2 in the DC approach. We see that, even 
though ohd vanishes linearly in x, it displays a much 
weaker x-dependence for x > 0.05. To make the last 
statement more quantitative consider the t/J — 3 re- 
sults, for which ohd{x — 0.05)/ max{a#£)(a;), < x < 
0.2} = 0.73 and a HD (x = 0.10)/max{o;rjjj(a:) ) < 
x < 0.2} = 0.93 in the SB framework, and olhd{x = 
0.05)/max{a_f/£i(a;), < x < 0.2} = 0.65 and otHD{x — 
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0.10)/max{a jF/D (a;),0 < x < 0.2} = 0.86 in the DC 
framework. We remark that the above behavior ap- 
pears to be a robust property of slave-particle formula- 
tions since it applies to two different slave-particle frame- 
works and various t/J values. This state of affairs should 
be contrasted with the approximate ps{0) oc x rela- 
tion applicable in both theories throughout the interval 
< x < 0.2 [Fig. 2(a)]. In addition to having different 
doping dependences, ps(0) and o-ud also display distinct 
parametric dependences on t/J in cither slave-particle 
approach. Hence, the way interactions and quantum fluc- 
tuations in doped Mott insulator superconductors renor- 
malize ps(0) differs from the way they renormalize a 13 . 
More importantly, we show that this difference is cap- 
tured by slave-particle approaches, which are often dis- 
missed on the grounds that they imply a ~ ps(0), a 
relation that counters experimental evidence 22-24 . Our 
calculation shows that this relation only holds in the 
asymptotic limit x — > 0, where the mismatch with exper- 
iments is expected since long-range AF order develops in 
material compounds. 

Fig. 2(b) further shows that in the hole doped cuprate 
regime, i.e. for t, t' , t" = tuD, t'liDi ^'hd^ a HD is approxi- 
mately a factor of 2.5 larger in the DC theory than in the 
SB approach [as expected from the quasiparticle current 
plots in Figs. 1(a) and 1(c)]. In particular, the inclu- 
sion of AF correlations, and the consequent momentum 
space anisotropy, brings cihd up to 0.5 < oihd J$ 0.6 
when 0.1 < x < 0.2 and t/J — 3, which is quantitatively 
consistent with experimental data 35-37 . 

The above AF correlations-driven enhancement of 
cxhd also renders thermally excited quasiparticlcs effec- 
tive in reducing the superfiuid stiffness in the DC the- 
ory. This effect is depicted in Fig. 2(c), which plots 

T Q C P = g T (t/J)T® p , where g T (t/J) = t/3J in the SB 
framework and grit/ J) = 2>J/t in the DC framework. 
In particular, for t,t',t" = tHD,t' HD ,t HD and x > 0.10, 
T^ p ~ J/10 in the DC theory, a value that is consis- 
tent with the cuprates' T c scale, and represents an order 
of magnitude improvement over the corresponding scale 
obtained within the SB approach. 

As Fig. 1(d) illustrates, AF correlations do not en- 
hance the nodal quasiparticle current in the electron 
doped regime and, in this case, the nodal current renor- 
malization factor ued is much smaller than its hole 
doped counterpart. This fact is attested by the small 
DC theory value of (cued / c^hd) 2 ~ 0.1 in Fig. 2(d), 
and is consistent with electron doped cuprates' exper- 
imental data showing a low temperature ps(T) hard 
to reconcile with gapless nodal excitations 38 ' 39 despite 
solid evidence for predominant rf x 2_ y 2-wave symmetric 
pairing 55-58 . This asymmetry between the electron and 
hole doped regimes is considerably larger in the DC 
approach than in the SB approach [Fig. 1(d)] and, 
therefore, our calculation suggests that the above appar- 
ent discrepancy between different experimental probes 
reflects the short-range AF correlations present in the 
strongly correlated SC state. 



The results in Fig. 2 disclose specific parametric de- 
pendences of the low temperature ps{T) on intermedi- 
ate values of t/J, which apply throughout a wide dop- 
ing range. These dependences can be remarkably dif- 
ferent in the SB and DC approaches, thus showing that 
they are strongly modified by the inclusion of local stag- 
gered spin correlations. For instance, the SB theory pre- 
dicts that anD decreases upon lowering t/J while the 
DC approach implies the opposite trend [Fig. 2(b)]. 
Only the DC theory result, however, correctly captures 
the well documented enhancement of quasiparticle fea- 
tures upon lowering t/J - it specifically predicts that 
«hj) oc ( J/t) 1 / 2 , a parametric dependence equal to that 
obtained for the nodal quasiparticle spectral weight in 
exact numerical calculations concerning the same param- 
eter regime 59 . Now consider the results in Fig. 2(c), 
which show that the quasiparticle-driven T c scale low- 
ers with increasing t/J in the SB approach, whereas it 
increases with t/J once AF correlations are included in 
the DC theory. The latter trend seems to be more con- 
sistent with experiments though. In fact, these support 
that the T c scale is, to a large extent, set by quasiparti- 
cles, and that superconductivity emerges in underdoped 
cuprates as a means to enhance the kinetic energy of 
charge carriers 60-62 (which is otherwise frustrated by the 
background staggered moment correlations). Hence, one 
expects T® p to grow with t/J, as obtained in the DC 
framework. 



IV. EXPERIMENTAL SIGNATURE OF AN 
APPLIED SUPERCURRENT 

Gauge invariance implies that, upon substituting A^ 
by the linear combination [A^ — (<f)j — <f>i)/2], where 
(f>i is the order parameter's phase, the Hamiltonians 
HfP F (Aij , aij) and H^ < F (A i j,a i j) describe the coupling 
between SC quasiparticles and an applied supercurrent, 
which can be addressed by experiments that probe single- 
electron physics. One such example is ARPES, which 
probes the single-electron energy dispersion E^^A), and 
that, at least in principle, provides the means to directly 
measure the quasiparticle current jb,k,u = ~ ■ 
Unfortunately though, such measurements require both 
good energy and momentum resolution, and are most 
likely unfeasible in underdoped cuprates, whose low en- 
ergy quasiparticle spectral features have small intensity 
and large widths. Alternatively, one may use STM, 
which has better energy resolution than ARPES. How- 
ever, STM is a local probe in real space and misses a con- 
siderable amount of momentum resolved information. In 
addition, since STM integrates over momentum space, it 
is only sensitive to the second power of an applied super- 
current's magnitude (as long as time-reversal symmetry 
remains unbroken). Still, as we show in what follows, 
STM can be used to probe certain qualitative features 
that derive from the underlying quasiparticle current mo- 
mentum space distribution. 
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A. Supercurrent dependence of tunneling 
conductance — BCS theory 

Wc first study the supercurrent dependence of the tun- 
neling conductance within the BCS theory. This allows 
us to introduce the general formal approach, as well as 
to estimate the (generic) order of magnitude of the effect 
produced by a supercurrent on the STM spectrum. 

The mean-field BCS superconducting Hamiltonian 

Hm C f S = ^2ci k e k c ak + ^2(c ak e af3 A k c,3,- k + h.c) (15) 



where e k is the normal state dispersion, A k is the gap 
function, and e a @ is the anti-symmetric tensor, can be 
rewritten as 

Hff F S = E a + J2(E% CS + e a k )bl k b ak , (16) 



where 



E 



BCS 



(4) 2 + |A fc p, 



(17) 



if we introduce the fermionic Nambu operators b ak 
Cffe = Ufe&Tfe + v kb\ _ h , c lk = u k b lk - v k b]- _ k (18) 

where u 2 k = |(1 + ■j&bs) and v k = |(1 - -slfer) ar e the 

BCS coherence factors. E in Eq. (16) is the ground- 
state energy, where the ground-state \SC) is determined 
by b ak \SC) = if E k + e a k > 0, and by bl k \SC) - if 



E k + e% < 0. 

From the above, we can obtain the density of states 
" lk T S(E k +el-e)ul 



N(e) 



b.z. (2tt) 
d 2 k 



B.Z. 



. (27T. 



2 S(E k + e% + e)v 2 . (19) 



In addition, since in layered materials like the cuprates 
we can approximately consider that the conductance be- 
tween a metal and a superconductor at a bias voltage V 
is proportional to the total tunneling density of states 
N(e) of the superconductor at e = V 63 ' 64 , we have 

d 2 k 



B.Z. 



. (2*) 



2 f(E k + e% + V)v 2 k (20) 



(i+e^) 2 and P = account for the effect 



where /(e) 

of thermal broadening 

We now apply the above formula to calculate the dif- 
ferential tunneling conductance between a metal and a 
simple BCS d-w&ve superconductor whose gap function 



k B T=10meV 




FIG. 3: A plot of [(#) k (V) - (#) Q (V)] / (#) (V 1/2 ) as 
a function of bias voltage V at three different temperatures 
and for na — 0.01. The solid line is for electron tunneling 
into the superconductor and the dashed line is for electron 
tunneling out of the superconductor. We use t — 0.3eV, /i = 
0.2eV, and A = 0.02eV which results in V gap ~ 0.066 volt. 
iw) (V) ~ iw)o(y} ^ as a temperature dependent peak 
around V = 0. The width of the peak at the half of maximum 
value is about 4fcsT. The height of the peak is proportional 
to l/T. 



is Afc = 2A[cos(fca;) — cos(k y )] , and whose dispersion 
in the absence of a supercurrent is e k = — 2i[cos(fc x ) + 
cos(fcy)] + /x. In the presence of a supercurrent J s in the 
plane, the energy dispersion shifts in momentum space 
as given by e k = c k+K , where k — A and where the vec- 
tor potential A is determined from London's equation 
J s = psA. Since the quasiparticle current is given by 
j k = — aE Q^ , the quasiparticle current distribution in 
the Brillouin zone controls the above energy dispersion 
shift. Consequently, the change in the tunneling spec- 
trum due to an applied supercurrent reflects the under- 
lying quasiparticle current distribution. 

To estimate the order of magnitude of the momentum 
shift At, we consider the effect of a supercurrent density 
of magnitude J s = 10 8 Amp/cm 2 on the surface of a su- 
perconductor. Such an current density can be achieved 
by passing 0.1 Amp of current through a superconducting 
thin film of 10/im wide and 0.01/im thick. The penetra- 
tion depth of the superconductor is Al = 0.1/im. (Note 
that, here, the London penetration depth is the pene- 
tration depth for a magnetic field perpendicular to the 
C u O plane.) Since the London penetration depth is 
given by Al = \J c 2 /Anpse 2 (the constants e and c are 
introduced for convenience), we have that 



4ttA 2 J s 



(21) 



where a embodies the effect of the quasiparticle current 
renormalization due to interactions. If we take the non- 
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interacting case, a = 1, and we find k = 1.75 x 10 5 /cm 
or na = 7 x 10~ 3 , where a = 3.8A is the lattice constant 
of the C u O plane. 

Let (jp 1 ) (V) be the differential tunneling conduc- 
tance in the presence of a supercurrent flowing in the 
C u O plane in the rr-direction. [(^) (V) then stands 
for the differential tunneling conductance in the ab- 
sence of a traversing supercurrent.] Fig. 3 plots 

[(w) K 00 ~ (#)o 00] I (#) (Vi/2) for the choice of 
parameters * = 0.3eV, fi = 0.2cV, A = 0.02eV, and 
na = 0.01. Here we use the scale factor (4^) (Vi/ 2 ), 
which denotes the differential tunneling conductance at 
Vi/2 = 0.033V = Vgap/2. Following the above calcula- 
tion we expect that, in an experimentally relevant con- 
text, the change in the tunneling spectrum is of the order 
of 0.001 — 0.01 of the original signal's magnitude. The 
effect of supercurrent on the tunneling dl/dV curve can 
also be studied by tunneling near a vortex. 



B. Supercurrent dependence of tunneling 
conductance — DC and SB theories 



We now focus on the particular case of a doped 
Mott insulator superconductor as described by the SB 
and DC theories. We thus extend previous calcu- 
lations of the tunneling differential conductance us- 
ing the SB and DC frameworks 44 ' 65 to account for 
the presence of an applied supercurrent. Specifically, 
we consider the mean-field Hamiltonians H^f F (Aij) = 
(A^aij)| 0= _( n sB)-i jjsb A and H MF (Aij) = 
(Aij , o>ij)\ a =-(ji DC )-\ ti dc A *° determine the de- 
pendence of the mean-field energy dispersions E^ B (A) 
and E^ c (A) on the gauge field A. We then straight- 
forwardly obtain the dependence of the differential tun- 
neling conductance (V) on A for both the SB and 
in the DC approaches as outlined in Sec. IV A. Let we, 
however, note a few technical differences between what 
follows and Sec. IV A. Firstly, below we assume that in- 
teractions provide the main contribution to the broaden- 
ing of spectral features of doped Mott insulators. Hence, 
instead of thermal broadening, we include the effect of 
a Lorentzian broadening parameterized by T = O.Oli (a 
value consistent with experiments 66 ' 67 ). Secondly, we ex- 
plicitly focus on the small K limit, in which case we can 
write 



h sb 

MF 
ttDC 
n MF 



dl_ 

dV 



(V) = 



dl 
dV 



(V)- 



n 



^)(^) + 0( K 4 ) (22) 



since in the absence of time-revesal symmetry breaking 
(w) K ( V ) = (w)- K (^)- Hence > below > we focus on the 
behavior of (M?) (V) instead of selecting a particular 
value of k. 

Figs. 4(a) and 4(b) depict the resulting DC and SB 
theory plots of 
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FIG. 4: (a) (j^)'L rm in DC theory. 
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bias gap 
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dl 



in SB the- 



ory, (c) (J£)'i in DC theory, (d) in SB theory. We 

plot the x — 0.05 (full line) and x = 0.15 (dashed line) results 
for t, t', t" — tHD,t' HD , t'x D , and use a Lorentzian broadening 
given by r = O.Oli. A finite T corresponds to a finite effetive 
temperature T whose value is propotional to and can be es- 
timated from the width of the peak at V = 0. The width at 
the half peak value is about 4fcsT (see Fig. 3). 



for x — 0.05 and x = 0.15 and in the subgap frequency 
range — V gap /2 < V < V ga p/2, where V is the bias voltage 
and V gap is the SC coherence peak voltage. (We only 
show results for the above values of V in order to focus on 
the specific experimental signature we discuss below.) In 
the above expression we normalize the second derivative 
of the tunneling conductance with respect to n so that it 
equals unity at V = 0. 

The interesting feature in Figs. 4(a) and 4(b) is that, 
out of all the curves in these two figures, the DC theory 
plot of ( ^ Ynorm f° r x = 0-05 stands out as the only curve 
which is clearly asymmetric around V = 0. To further 
emphasize the asymmetry in the DC theory x = 0.05 
curve, as well as the symmetry around V — of all other 
curves, in Figs. 4(c) and 4(d) we plot (4k)'L, where 



dl\" , s ( dl\' s 
dv) { - v) -\dv) in 



(24) 



The question then arises of what the physical reason 
that singles out the DC theory x — 0.05 curve is. We find 
there are at least three reasons to associate the tilting to- 
ward the negative bias side in the DC theory x = 0.05 
curve to the presence of local staggered moment corre- 
lations. Firstly, such a qualitative feature is altogether 
absent in the SB results. Secondly, the above asymme- 
try develops upon lowering x, which is known to enhance 
the signatures of AF correlations. Lastly, it naturally fol- 
lows from the DC theory two-band picture that describes 
the interplay between coexisting AF and SC correlations 
at short length scales 31 ' 32 . To clarify the latter point, 
we remark that the DC mean-field theory contains two 
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different families of fermions, namely spinons and do- 
pons, whose dispersions are determined by Eqs. (9) and 
(10). Applying a supercurrent shifts the spinon and do- 
pon bands relatively to each other and, thus, affects the 
electronic spectral weight transfer to low energy [which is 
determined by the hybridization of spinons and dopons 
described in Eq. (11)]. This spectral weight transfer 
is reduced mainly in those regions of momentum space 
where the second derivative with respect to momentum 
of the energy difference between both bands is larger, 
which happens to occur close to the peak of the AF-likc 
dopon dispersion, hence close to (7r/2,7r/2). Since the 
spinon nodal point shifts away from (7r/2,7r/2) toward 
(0,0), the above spectral weight reduction is stronger in 
the positive bias side, as obtained in Fig. 4(a) [this argu- 
ment implies that if the nodal point were to shift toward 
(71", 7r), the curve would rather tilt in the opposite 

direction] . 

From the above argument, the DC theory bias asym- 
metry in (4^)" relies on two things, namely, the nodal 
point shift away from (7r/2,7r/2) and the presence of 
strong local AF correlations. There exists ample exper- 
imental evidence for the former 68-71 . As to the latter, 
AF correlations were proposed to underlie the momen- 
tum space anisotropy that weakens the differential tun- 
neling conductance SC coherence peaks 44 . Therefore, 
we propose that if, indeed, local AF correlations are the 
cause of the aforementioned differentiation of the nodal 
and antinodal momentum space regions, an asymmetry 
around V = should be detected in the (4^)" curve mea- 
sured by STM experiments in the large gap (and small 
coherence peak) regions of inhomogeneous BSCCO sam- 
ples traversed by a supercurrent. This effect should be 
weaker, if at all observable, in the small gap (and large 
coherence peak) regions of these same samples. 

We remark that the above asymmetry in is max- 

imal at low values of V (in the above calculation, this 
asymmetry peaks around V ~ 0.15V gap )- At this en- 
ergy scale the differential conductance in the absence of 
a supercurrent, namely (V), is nearly symmetric 
around V = 0, which should facilitate the detection of 
the aforementioned low energy bias asymmetry. The low 
bias tunneling spectrum is also nearly spatially homoge- 
neous, a fact that should also facilitate the detection of 
the spatially inhomogeneous asymmetry of (JpO" a * l° w 
bias. In this context, we find that, within the DC frame- 
work, the low bias (7^7)" spatial inhomogeneity correlates 
with three aspects of the differential tunneling conduc- 
tance at higher energy, namely, the gap size, the size of 
SC coherence peaks, and the high energy asymmetry be- 
tween positive and negative bias. 

The above proposal of a specific experimental signa- 
ture of a traversing supercurrent in the STM spectra of 
BSCCO samples relies on a calculation that assumes an 
homogeneous system. There are two reasons to believe 
our proposal is robust to the BSCCO samples' spatial 
inhomogeneity. Firstly, the large BSCCO 's STM spec- 
tral diversity can be reproduced in homogeneous systems 
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proximate to a Mott insulator transition 44 . Secondly, the 
spatial inhomogeneity correlates with off-plane disorder 72 
which affects the effective parameters t' and t" , but not 
t and J 73 . In Fig. 5 we depict the DC and SB theory 
results for the curves (4^)- with x = 0.05 and x = 0.15, 



-2t" = 0M' HD and t' 



-It" = 2t', 



HD- 



as well as t' 

These show that our (4^)'L results are almost insensitive 
to changes in t' and t" in the range 2t' HD < t! < 0.5t' HD 
and Q.bt'ljjj < t" < 2t" HD , which argues in favor of the 
robustness of the experimental effects discussed above. 



V. SUMMARY 

As discussed in the introduction, experimental data 
together with theoretical arguments support the impor- 
tant role of thermally excited SC quasiparticles in set- 
ting the T c scale of underdoped cuprates. In this paper, 
we use two different wave functions, namely the SB 27,28 
and DC 31,32 d-w&ve SC wave functions, together with 
the tt't" J model Hamiltonian, to show that the com- 
bined effect of slave- particles and local AF correlations 
reproduces non-trivial aspects of these quasiparticles' low 
energy and long wavelength electromagnetic response. 

Slave-particle formulations are attractive in that they 
provide a microscopic description of doped Mott insula- 
tor superconductors which yields both (i) a non- vanishing 
quasiparticle d-wave gap and (ii) a vanishing effective 
density of charge carriers as the half-filling composition 
is approached. However, previous work on the SC electro- 
magnetic response in slave-particle frameworks 9 was not 
consistent with a third crucial experimental fact, namely, 
that the nodal current renormalization factor a displays 
a much weaker x-dependence than ps (0) . The aforemen- 
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tioncd work was concerned with the x — > limit, where 
real materials display long-range AF order and do not su- 
perconduct. Even though variational Monte-Carlo stud- 
ies have extended the calculation of a to doping values 
beyond the above limit 74 , this technique only obtains a 
maximum bound on the value of ps(0) 75 . In this paper, 
we relax the no-double occupancy constraint (which is 
implemented exactly in the variational Monte-Carlo ap- 
proach) and only include (the gapped) gauge fluctuations 
at the Gaussian level. This allows us to calculate the 
temperature dependent electromagnetic response in the 
static and uniform limit for two different slave-particle 
approaches and, consequently, we are able to compare 
the doping dependence and the t/J parametric depen- 
dence of both ps(0) and a. Interestingly, we find that, 
even though a ~ ps(0) ~ x in the limit x — > 0, away 
from this limit a is much more weakly rr-dcpendcnt than 
ps(0). This result applies to both utilized slave-particle 
frameworks and, we propose, may be generic to slave- 
particle formulations. 

In this paper, we also compare the SB and the DC 
theory results to learn how high energy and short-range 
staggered moment correlations affect the SC electromag- 
netic response in the static and uniform limit. We find 
that inclusion of these correlations improves SB results, 
as we summarize below. For instance, the SB and the 
DC theories imply different parametric dependences on 
intermediate t/J values and, in Sec. IIIB, we argue in 
favor of the DC theory expectations. We also find that 
local AF correlations, which enhance the quasiparticle 
current in specific momentum space regions that differ 
for the hole and electron doped regimes, provide a micro- 
scopic rationale for the experimentally observed hole vs. 
electron doped asymmetry in the nodal quasiparticlcs' 
electromagnetic response. The aforementioned quasipar- 
ticle current renormalization further brings quantitative 



agreement with hole doped cuprates' experimental data, 
specifically a and the T c scale for x > 0.10, if we use 
physically relevant bare parameters in the tt't" J model 
Hamiltonian. To understand the above improvement 
upon inclusion of AF correlations note that, in the SB 
framework, the projection constraint is implemented af- 
ter integrating out the 5(7(2) gauge field, a procedure 
that enhances this type of AF correlations 76 . To check 
the consistency of this interpretation, we refer to a varia- 
tional Monte-Carlo study 74 that exactly enforces the no- 
double occupancy constraint on BCS wave functions and 
whose values of a are indeed larger than those we obtain 
in this paper's SB approach. 

Finally, we point out that a momentum dependent, and 
thus energy dependent, coupling of quasiparticles to an 
electromagnetic gauge field is reflected in the local den- 
sity of states in the presence of an applied supercurrent. 
In this context, we derive (DC theory specific) qualitative 
predictions for the effect of short-range staggered corre- 
lations between local moments in the tunneling spectrum 
of a sample traversed by a supercurrent, namely, this the- 
ory implies an asymmetry in ( j£)"(V) around V = for 
underdoped samples. This effect may be probed by STM 
experiments, which thus could distinguish the DC theory 
description of the high-T c superconductors from the SB 
theory and the conventional BCS theory. 
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